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ABSTRACT 


Various  computational  methods  and  operational  computer  codes  used  to  predict  the 
aerodynamic  coefficients  and  separation  trajectories  of  aircraft  stores  are  examined.  The 
semi-empirical  aeroprediction  code  Missile  DATCOM  is  used  to  obtain  the  coefficients 
of  a  modeled  store.  These  coefficients,  together  with  the  modeled  ejection  forces,  are 
used  in  free-stream  state-space  equations  of  motion  to  predict  the  store  trajectory.  The 
results  are  compared  with  the  Nielson  Engineering  and  Research  (NEAR)  store  separation 
code  which  provides  accurate  trajectory  profiles,  for  speeds  below  the  critical  speed,  by 
use  of  a  vortex-lattice  and  panel  method.  Modification  of  the  Missile  DATCOM 
aerodynamic  coefficients  provides  single-point  state-space  prediction  of  the  store  pitch 
trajectory  within  30%  of  the  NEAR  code  results.  Store  trajectories  were  restricted  to  the 
first  0.2  seconds  of  free  flight. 
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I .     STORE  SEPARATION  INTRODUCTION 

A .   BACKGROUND 

The  prediction  of  the  trajectory  of  a  store  ejected  from 
an  aircraft  is  of  major  concern  to  defense  aviation.  The 
increasing  requirement  for  aircraft  to  fulfill  a  multiple  role 
in  the  hostile  environment  demands  that  it  also  carry  an  ever- 
increasing  variety  of  stores.  These  stores  range  from 
missiles  and  bombs  stockpiled  for  many  years  to  new, 
aerodynamically  complex  stores  with  lifting  bodies. 

The  aim  of  this  work  is  to  review  current  methods  of 
predicting  store  releases,  compare  results  of  two  different 
methods,  and  to  present  an  easy-to-use  methodology  for 
investigating  the  subsonic  release  of  a  unsophisticated  store. 

A  non-axisymmetric  pod  ejected  from  the  F/A-18  outboard 
pylon  is  presented  as  an  example.  This  store  was  chosen 
because  it  is  of  practical  value  due  to  the  extensive  use  of 
these  modified  pods  at  the  Pacific  Missile  Test  Center,  at 
Point  Mugu,  California.  Non-axisymmetrical  stores  are  very 
difficult  to  model  and  present  a  special  problem  in  the  store 
separation  field.  This  pod  is,  therefore,  modeled  as 
axisymmetric  and  non-axisymmetric,  and  the  results  are 
compared.  In  addition,  the  store  lacks  aerodynamic  control 
surfaces   and   is   highly   unstable.     Application   of   the 


methodology  presented  here  should  eventually  provide  the  user 
with  accurate  results  for  a  minimal  time  and  cost  expenditure. 

B.   MATHEMATICAL  MODELING 

There  are  many  ways  to  mathematically  model  store 
separation.  Economics  and  accuracy  are  always  the  main 
considerations  when  exploring  the  codes  and  methods  to  be 
used.  For  simple  cases,  intermediate  approaches  such  as  panel 
methods  and  solutions  of  the  Euler  equations  provide 
sufficient  accuracy.  Calculations  involving  flow  separation 
and  shock  wave  interference  become  too  complex  for 
intermediate  methods  and  acceptable  results  can  only  be 
obtained  using  modern  computational  fluid  dynamic  (CFD) 
techniques  involving  solutions  of  the  Navier-Stokes  equations 
[Ref.  l:p.  7-2].  For  speeds  below  the  subsonic  Mach  critical 
speed,  component  buildup  methods  (empirical  and  semi- 
empirical)  ,  and  panel  methods  provide  the  necessary  accuracy 
and  are  much  easier  to  use.  The  codes  chosen  for  this  study 
are  explained  in  detail  in  their  respective  sections . 

The  F/A-18  model  was  obtained  from  Mr.  L.  L.  Gleason  of 
the  Ordnance  Systems  Department,  Naval  Weapons  Center,  at 
China  Lake,  California.  The  F/A-18  was  chosen  due  to  its  wide 
variety  of  store  carriage  and  also  because  of  its  solid  future 
with  Naval  aviation. 

The  store  model  is  based  upon  both  the  Missile  Datcom 
method  [Ref  2]  and  also  by  the  Nielson  program  [Ref  3]  .   These 


will  be  discussed  in  Sections  III  and  IV.  Appendix  A  contains 
the  physical  description  of  the  pod  and  its  inertial 
characteristics . 

C .   METHODOLOGY 

The  initial  investigation  (Section  II)  consisted  of 
looking  at  the  ejector  forces,  combined  with  the  store's 
inertial  characteristics,  to  give  a  sense  of  the  magnitude  of 
initial  velocities  and  moments.  These  forces  and  moments, 
while  easy  to  predict,  were  insufficient  to  provide  real 
evidence  of  a  safe  release  or  of  a  release  problem.  A  large 
pitch-down  moment  or  a  large  vertical  velocity  may  be  off-set 
by  unforeseen  aerodynamic  forces .  The  velocities  and  moments 
calculated  will  be  used  in  Section  III  and  compared  with 
results  of  Section  IV  for  trajectory  prediction. 

Section  III  introduces  a  method  of  predicting  the  store's 
longitudinal  and  lateral  aerodynamic  coefficients  using  a 
missile  code  (Missile  DATCOM)  developed  by  McDonnell  Douglas 
Missile  Systems  Company  for  the  Flight  Dynamics  Laboratory 
(FDL) ,  Wright-Patterson  Air  Force  Base  (WFAFB) ,  Ohio  [Ref .  2]  . 
These  aerodynamic  coefficients  are  then  used  in  longitudinal 
and  lateral  equations  of  motion  to  predict  the  free-stream 
trajectory  of  the  store.  A  brief  explanation  of  the 
derivation  of  these  equations  of  motion  is  provided  to  the 
reader.  An  accurate  representation  of  the  aerodynamic 
characteristics  is  essential  for  the  full  prediction  of  the 


shape's  trajectory.  Missile  DATCOM  is  based  on  the  body 
buildup  method  and  includes  a  number  of  prediction  methods  for 
each  component  of  the  configuration.  Other  codes,  such  as 
MISL3,  are  also  available,  each  with  its  own  relative 
strengths . 

Section  IV  describes  a  computer  prediction  method 
developed  by  Nielson  Engineering  and  Research,  Inc.,  (NEAR), 
also  under  contract  to  the  United  States  Air  Force.  This  code 
provides  a  six  degree-of-f reedom  (6D0F)  simulation  which  takes 
into  account  the  aircraft  flowfield  using  vortex-lattice  and 
panel  methods.  The  NEAR  simulation  provides  a  high  fidelity 
representation  for  the  subsonic  case,  but  also  requires  the 
most  effort  and  computer  capability.  An  attempt  has  been  made 
to  minimize  the  complexity  of  the  input  to  the  program.  Once 
a  good  data  base  has  been  developed  the  NEAR  code  should  prove 
easy  to  use.  It  has  been  used  extensively  throughout  the 
defense  industry,  and  has  undergone  many  modifications  to 
incorporate  improvements  and  options. 

Section   V   compares   the   trajectory   results   of   the 
linearized  aerodynamic  simulation  method  and  the  NEAR  code. 
A  modification  to  the  linear  coefficients  used  in  the  linear 
aerodynamic  simulation  is  also  discussed. 


II.   EJECTION  FORCES  AND  MOMENTS 

Although  quite  uncomplicated  in  nature,  calculating  the 
ejection  forces  and  moments  provides  the  user  with  a  feel  for 
the  magnitude  of  the  initial  movement  of  the  store.  This 
procedure  should  be  done  prior  to  using  any  other  method  as  a 
preliminary  step  in  order  to  prepare  for  the  simulation. 

A.   EJECTOR  CHARACTERISTICS 

The  non-aerodynamic  forces  involved  in  the  ejection  of  a 
store  depend  upon  the  ejector  cartridge,  the  ejector  rack,  the 
weight  of  the  store,  and  the  rigidity  of  the  wing.  Each  of 
these  parameters  determine  the  resultant  moments  and  forces 
provided  to  the  store.  The  store  and  the  wing  (or  fuselage) 
were  considered  as  rigid,  thus  neglecting  any  aeroelastic  or 
structural  bending  effects.  None  of  the  methods  discussed 
here  address  aeroelastic  bending  due  to  the  added  complexity 
of  the  problem. 

The  Douglas  production  BRU-32/A  Bomb  Ejector  Rack  combines 
two  sets  of  hooks,  one  set  at  14-inch  spacing  and  one  set  at 
30-inch  spacing,  with  an  ejection  system  designed  for  carriage 
of  stores  with  suspension  lugs  per  MIL-A-8591.  Two 
electrically  initiated  CCU-45B  cartridges  are  used  for  store 
release  and  ejection.   The  self-retracting  ejector  pistons, 


spaced  symmetrically  20  inches  apart  at  each  end  of  the  rack, 
have  a  piston  stroke  of  six  inches.  These  pistons  are  spring 
loaded  against  the  store  during  loading  to  prevent  impact  of 
pistons  during  firing.  The  orifice  sizes  can  be  varied,  by 
replacement,  to  provide  force  and  pitch  control  for  store 
separation.  [Ref.  4] 

The  best  source  of  ejector  data  is  from  the  Aircraft 
Ordnance  Procedures  (AOP)  contained  in  the  Aircraft  Stores 
Interface  Manual,  (Reference  4)  .  This  manual  contains  a 
complete  description  of  the  rack  under  consideration.  Here, 
the  example  separation  uses  two  14-inch  spaced  hooks  and 
0 . 118-Diameter  orifices.  The  force  diagram  corresponding  to 
this  orifice,  modeled  from  the  AOP,  is  shown  in  Figure  1. 
Sufficient  data  points  were  read  from  the  AOP  graph  in  order 
to  curve-fit  the  points  using  Computer  Associates' 
CricketGraph*  software  [Ref.  5].  This  provided  two  things; 
1 .  )  integration  of  the  polynomial  equation  gave  the  total 
impulse  value,  and  2.)  Nielson  software  in  Section  IV  uses  a 
fifth-order  polynomial  for  calculations.  Note  that  the  AOP 
shows  the  total  force  provided  by  the  ejector,  while  the  NEAR 
simulation  requires  a  force-per-e jector-foot  polynomial. 
Given  the  shape  of  the  force  curve,  a  triangular 
representation  of  the  curve  gave  a  value  of  348  lbf-sec  versus 
the  integrated  value  of  344.9  lbf-sec,  which  corresponds  to  an 
error  of  only  1%.  Therefore  elementary  estimates  of  the  total 
impulse  are  practical. 
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Figure  1    Ejector  Force 


B .   FORCES 

1 .  Theory 

Using  the  conservation  of  linear  momentum  theory 
presented  in  Appendix  B,  the  pod  was  modeled  as  a  simple  beam. 
The  weight  was  concentrated  at  the  center  of  gravity  and  the 
resultant  ejection  forces  were  placed  at  the  ejector  feet 
locations.  The  ejection  force  time  histories  are  treated  as 
an  equivalent  force  impulse  which  results  in  a  change  in 
linear  momentum  as  shown  by  the  following  equation: 
JlFz  dt  =  G2  -  G1   =  (M*VS)  2  -  (M*V.)  J 

The  integral  for  representing  the  linear  impulse  is 
based  upon  the  fifth-order  polynomial  force  time  history  shown 
here . 

y=  280. 2-608022. 8x  +  180688545 . 8x2  -  10638451765 . 5x3  + 
234526098234. 3x4  -  1783873760363 . 4x5 

2 .  Example 

For  our  example,  the  BRU-32A  bomb  rack  used  two  CCU- 
45B  cartridge-activated  devices  (CADS) .  The  peak  force  was 
14,500  lbf  and  the  pulse  duration  was  50  msec.  The  total 
linear  impulse  calculated  was  345  lbf-sec .  Using  a  weight  of 
371  pounds  for  the  pod,  the  end-of-stroke  velocity  was  29.9 
feet  per  sec  (fps)  .  This  value  was  used  in  Section  III  as  the 
initial  velocity  for  the  trajectory  simulation. 


The  graphs  contained  in  the  AOP  shows  end-of-stroke 
velocity  for  a  371  lb  store  to  be  approximately  24.5  f  ps . 
This  graph  is  derived  from  empirical  data  and  should  provide 
the  most  accurate  representation  of  the  actual  release 
velocity.  The  values  of  24.5  and  29.9  fps  are  compared  with 
the  end-of-stroke  store  velocity  of  30.9  fps  predicted  from 
the  NEAR  code  in  Section  IV. 

C .   MOMENTS 

1 .  Theory 

The  pod  was  modeled  as  a  simple  beam.  The  inert ial 
characteristics  of  the  store  are  listed  in  Appendix  A.  The 
two  ejector  feet  provide  the  beam  with  a  moment.  The 
equations  used  are  elementary,  and  with  certain  assumptions 
the  angular  velocity  can  be  calculated.  The  equations  are 
based  on  conservation  of  angular  momentum.  The  equations  and 
assumptions  are  presented  in  Appendix  B  for  completeness.  The 
resulting  equation  is  shown  here. 

jlMy  dt  =  iyy  *  (^ 

2 .  Example 

The  center  of  gravity  was  offset  aft  of  the  center  of 
the  ejector  feet  by  4.05  inches.  The  impulse  thus  provided 
the  pod  with  1.12  radians  per  second  (rps)  or  64.1  deg/sec 
initial  angular  velocity.  The  question  now  is,  are  these 
reasonable  values?  At  the  end  of  0.429  seconds,  the  shape 
will  be  one  body-length  below  the  aircraft.  At  the  same  time, 


the  shape  will  have  rotated  27.5  degrees  nosedown.  Therefore, 
without  considering  the  aerodynamics  of  the  vehicle,  the  pod 
seems  to  have  cleared  the  wing.  A  good  rule  of  thumb  is  2-3 
body  lengths  clearance,  although  exceptions  do  occur. 


D.   PRELIMINARY  OBSERVATION 

This  initial  calculation  provides  an  approximate  estimate 
of  the  forces  and  moments  involved.  The  values  obtained  are 
of  sufficient  accuracy  for  calculations  such  as  end-of-stroke 
loading,  etc.  Excessively  large  separation  velocities  are 
usually  an  indication  of  a  miscalculation.  Translational 
velocities  should  range  from  10-30  fps  and  pitching  velocities 
from  0-6  rps . 
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III.   AERODYNAMIC  CALCULATION 

A .   THEORY 

Many  physical  systems  can  be  modeled  by  second-order 
differential  equations.  The  mathematical  treatment  of  fixed- 
wing  flight  vehicle  motions  was  first  developed  by  G.H.  Bryan. 
He  laid  the  mathematical  foundation  for  airplane  dynamic 
stability  analysis,  developed  the  concept  of  the  aerodynamic 
stability  derivative,  and  recognized  that  the  equations  of 
motion  could  be  separated  into  a  symmetric  longitudinal  motion 
and  an  unsymmetric  lateral  motion.  Experimental  studies  were 
initiated  by  L.  Bairstow  and  B.M.  Jones  of  the  National 
Physical  Laboratory  in  England,  and  by  Jerome  Hunsaker  of  the 
Massachusetts  Institute  of  Technology  to  determine  estimates 
of  the  aerodynamic  stability  derivatives  in  Bryan's  theory. 
In  addition  to  determining  stability  derivatives  from  wind- 
tunnel  tests  of  scale  models,  Bairstow  and  Jones 
nondimensionalized  the  equations  of  motion  and  showed  that, 
with  certain  assumptions,  there  were  two  independent 
solutions,  i.e.,  one  longitudinal  and  one  lateral.  These  two 
solutions  provide  the  free-stream  trajectories  we  seek.  [Ref . 
6:p.  113] 

The  dynamic  stability  characteristics  of  a  pod  can  be 
represented  by  six  equations  of  motion,  three  for  the  forces 
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involved  X,  Y,  and  Z  and  three  for  the  moments  L,  M  and  N. 
The  force  equations  relate  the  forces  acting  on  the  body  to 
the  corresponding  linear  accelerations  and  the  moment 
equations  relate  the  moments  to  the  corresponding  angular 
accelerations.  It  is  usually  possible  to  consider  the 
longitudinal  motions  completely  separately  from  the  lateral- 
directional  motion,  by  neglecting  the  various  coupling  terms. 
[Ref.  7:p.  14] 

As  a  caveat  to  the  use  of  this  method,  these  equations  are 
the  linearized  version  and  are  only  valid  up  to  approximately 
10  degrees  angle-of-attack .  Treatment  of  the  non-linear 
aerodynamic  coefficients,  while  not  extremely  difficult,  does 
require  knowledge  of  the  behavior  of  the  coefficients.  Since 
the  goal  of  this  report  is  to  predict  the  behavior  of  new 
shape  configurations,  this  knowledge  is  not  presumed. 
Therefore,  the  results  obtained  in  this  fashion  are  only  valid 
for  estimating  the  motion  in  the  first  fractions  of  a  second. 
These  results  are  compared  with  those  obtained  in  the  more 
rigorous  method  of  Section  IV. 

1 .  Longitudinal  Equations 

The  rigid-body,  longitudinal  equations  of  motion  can 
be  developed  from  Newton's  second  law: 

^Pitching  moment 3=2^^=1^*0 
The  assumption  that  body  motion  consists  of  small  deviations 
from  its  equilibrium  flight  condition  allows  us  to  use 
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perturbation  theory  to  examine  the  aerodynamic  force 
derivatives  in  terms  of  angle-of-attack  (AOA) ,  vertical 
velocity,  pitch  angle,  and  pitch  rate,  by  means  of  a  Taylor 
series  expansion.  The  X-force,  Z-force,  and  pitching  moment 
equations  comprise  the  longitudinal  equations .  To  separate 
these  equations  from  the  lateral  equations,  they  must  not  be 
coupled.  This  is  a  reasonable  assumption  given  the  nominal 
geometric  shape  of  missiles  or  pods  [Ref .  6]  .  The  body  is 
constrained  to  move  in  a  vertical  plane  and  is  free  to  pitch 
about  its  center  of  gravity.  For  a  comprehensive  derivation 
of  the  equations  of  motion,  see  Reference  6.  The  resulting 
equations  of  motion  are  shown  in  Appendix  C. 
2 .  Lateral  Equations 

The  Y-force,  rolling,  and  yawing  moment  equations 
comprise  the  lateral  equations  of  motion.  Once  again,  these 
equations  are  uncoupled  from  the  longitudinal  equations  by  the 
assumption  of  small  cross-component  moments  of  inertia.  That 
is,  Ixy,  Ixz,  Iyx  are  small  in  comparison  to  the  principal  axis 
moments.  These  equations  are  also  valid  assuming  only  small 
variations  in  displacement  and  velocity.  The  lateral 
equations  are  derived  from  the  following  Newton's  laws: 

ZRolling  moment s= 1^*0 

ZYawing  moment s= I „*\|/ 
The  third  equation  is  derived  from  the  Taylor  series  expansion 
of  the  side  force  derivative.   The  lateral  motion  of  the  pod 


13 


disturbed   from   its   equilibrium   state   is   a   complicated 
combination  of  rolling,   yawing,   and  sideslipping  motions. 
Assuming  the  cross  products  of  inertia  are  ignored,  some  of 
the  coupling  terms  can  be  simplified. 
3  .  State  Variable  Representation 

The  linearized  longitudinal  and  lateral  equations 
developed  above  are  simple,  ordinary  linear  differential 
equations  with  constant  coefficients.  The  coefficients  in  the 
differential  equations  are  made  up  of  the  aerodynamic 
stability  derivatives,  mass,  and  inertia  characteristics  of 
the  pod.  These  equations  can  be  written  as  a  set  of  first- 
order  differential  equations  in  the  state-space  (state 
variable)  form: 

{x}  =  [A]*{x}  +  [B]*{T1} 

where  {x}  is  the  state  vector,  {T|}  is  the  control  vector  and 
the  matrices  [A]  and  [B]  contain  the  pod's  dimensional 
stability  derivatives.  In  our  case  the  pod  does  not  have  any 
active  control  surfaces.  For  missile  launches,  however,  the 
control  vector  and  the  B-matrix  would  be  used  to  represent  the 
control  surfaces  and  control  input.  The  simplified  state- 
space  representations  of  the  longitudinal  and  lateral 
equations  of  motion  are  shown  in  Appendix  C.   The  required 
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coefficients  listed  in  Appendix  C  are  derived  using  the 
procedures  in  Section  III.B. 

The  state-space  representation  of  the  equations  of 
motion  can  be  solved  simultaneously  using  matrix  software  such 
as  MATLAB*  [Ref .  8]  or  Control-C*  [Ref .  9]  .  Given  the 
conditions  of  the  ejection,  namely,  initial  translational 
velocities,  initial  angular  velocities,  and  attitude  angles, 
solutions  of  the  state-space  equations  can  be  used  to  predict 
the  trajectory  of  the  pod.  The  initial  angles  and  velocities 
are  the  initial  conditions  imposed  upon  the  state-space 
equations.  The  Control-C*  commands  for  this  procedure  are 
contained  in  Appendix  D . 
B.   MISDATCOM 

MISDATCOM  was  developed  by  McDonnell  Douglas  Astronautics 
Company,  St  Louis,  Missouri,  for  the  Flight  Dynamics 
Laboratory  of  the  Air  Force  Wright  Aeronautical  Laboratories, 
Wright -Patterson  AFB .  The  program  was  completed  in  December 
1985.  The  Missile  Datcom  was  created  to  provide  an 
aerodynamic  design  tool  which  has  the  predictive  accuracy 
suitable  for  preliminary  design,  and  the  capability  for  the 
user  to  easily  substitute  methods  to  fit  specific 
applications.  [Ref.  2:p.  1] 

1 .  Theory 

There  are  many  different  methods  to  predict  the 
aerodynamic  static,  dynamic,  and  control  characteristics  of 
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missiles.  Component  build-up  was  chosen  as  the  most  suitable 
for  this  program.  Although  panel  methods  are  better  suited 
for  arbitrary  configurations,  component  build-up  was  chosen 
due  to  the  accuracy  provided  for  conventional  configurations 
and  for  the  ability  to  do  parametric  studies  easier. 

Basically,  component  build-up  consists  of  using 
various  methods  to  compute  the  characteristics  (skin  friction, 
force  and  moment  coefficients,  panel  loading,  .  .  .)  of  the 
individual  configuration  components.  The  various  methods  are 
chosen  for  their  applicability  to  the  configuration  or  flight 
condition.  Then  the  components  are  combined.  Previous 
methods  of  combining  the  components  (fins,  body,  engine  inlet) 
consisted  of  adding  the  individual  coefficients  and  then 
multiplying  the  sum  by  some  interference  factor  obtained  using 
slender  body  theory.  The  approach  t~ken  with  MISDATCOM  was  to 
use  the  "equivalent  angle  of  attack  method"  developed  by 
Nielson  Engineering  and  Research,  Inc.  (NEAR) .  This  method 
assumes  that  the  panel  loading  for  a  given  panel  angle-of- 
attack  is  unique.  With  this  method  the  panel  angle  of  attack 
is  computed  including  the  effect  of  panel  roll  orientation 
with  respect  to  the  free  stream  velocity  vector,  panel 
proximity  to  the  fuselage  or  to  other  panels,  and  external 
vortex  flow  field  effects.  Then  the  isolated  panel 
characteristics  are  interpolated  at  the  panel  equivalent  angle 
of  attack  to  yield  the  panel  load  when  mounted  on  a  body  in 
combination  with  other  surfaces.  [Ref.  2] 
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2 .  Procedure 

The  procedure  for  using  MISDATCOM  is  straightforward. 
The  pod  shape  is  modeled  using  simple  geometric  shapes.  The 
previous  release  of  MISDATCOM  software  could  only  handle 
axisymmetric  shapes  or  forms  with  a  vertical  plane  of 
symmetry.  Unfortunately,  the  pod  under  consideration  is  non- 
axisymmetric  and  therefore  the  results  are  not  entirely  valid. 
The  latest  version,  however,  can  model  some  non-axisymmetry 
through  the  use  of  a  "protuberance"  option.  This  is  the  April 
1991  release  and  is  now  available. 

For  this  investigation,  a  simulation  run  was  made 
using  both  versions  and  a  comparison  of  the  aerodynamic 
coefficients  was  made.  The  pod  radius  was  also  varied  from 
the  minimum  pod  to  the  maximum,  and  the  resulting  coefficients 
compared.  The  MISDATCOM  code  coefficients  were  somewhat 
insensitive  to  radius  changes  of  this  magnitude,  therefore, 
the  minimum-radius,  axisymmetric  case  was  retained  for 
comparison  with  the  minimum-radius  non-axisymmetric  case. 
Figure  2  shows  the  pod  semi-profile  of  the  MISDATCOM  procedure 
alongside  the  pod  semi-profile  of  the  NEAR  simulation.  Figure 
3  contains  a  sample  input  to  the  program  and  Figure  4  lists  a 
partial  sample  output.  The  actual  input  and  output  of  the 
program  are  contained  in  Appendix  E.  These  aerodynamic 
coefficients  are  next  used  in  the  longitudinal  and  lateral 
state-space  equations  of  motion  to  predict  the  store 
trajectory . 
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Figure  2    MISDATCOM  Pod  Shape 
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Figure    4         MISDATCOM   Output 
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C.   AERODYNAMIC  CALCULATION 
1 .  Longitudinal 

The  aerodynamic  coefficients  derived  from  the 
MISDATCOM  were  entered  into  the  plant  matrix  [A]  of  the 
longitudinal  equations  of  motion.  The  initial  conditions  were 
applied  and  time  history  of  the  pitch  angle,  pitch  angular 
velocity,  and  vertical  velocity  was  found.  Figures  5-8  show 
that  the  pod  assumes  a  large  nose-down,  highly  divergent 
pitching  motion.  This  is  understandable  since  there  are  no 
control  surfaces  to  make  the  pod  stable.  Examination  of  the 
roots  of  the  plant  matrix  [A]  indicates  an  unstable  flight 
vehicle.  The  pitch  angle  and  pitch  rate  shown  are  in  the 
store  body  coordinate  frame. 

These  results  are  valid  only  for  the  range  of  linear 
values  of  the  aerodynamic  coefficients.  This  is  approximately 
up  to  10-12  degrees  angle-of-attack.  Therefore,  an  estimate 
of  the  non-linear  behavior  of  the  pod  after  an  AOA  of  10 
degrees  is  reached  is  necessary.  Due  to  the  large  pitch 
angular  velocity,  it  is  obvious  that  this  limit  is  reached 
after  only  0.2  seconds.  Section  V  discusses  the  acceptable 
time  range  where  these  results  are  valid. 

To  estimate  the  non-linear  behavior,  the  dominating 
terms  in  the  [A] -matrix  must  be  determined.  This  provides  an 
insight  into  the  possible  range  of  values  to  substitute  into 
the  [A] -matrix.   A  discussion  of  a  possible  approach  to  this 
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problem  is  presented  in  Section  V.  To  completely  investigate 
this  area  is  beyond  the  scope  of  the  method  outlined  here. 
Comparison  of  the  "linear"  trajectory  obtained  here  will  be 
made  with  the  trajectory  obtained  in  Section  IV,  which  does 
take  into  account  the  non-linearity  of  the  coefficients. 
2 .  Lateral 

Due  to  the  nearly  vertical  forces  and  moments  provided 
by  the  ejection  rack  during  straight  and  level  flight,  the 
free-stream  lateral  equations  of  motion  will  not  provide  us 
with  any  insight  into  the  safe  jettison  of  the  shape. 
However,  if  the  pod  were  experiencing  sideslip,  then  the 
sideslip  could  be  entered  into  the  state-space  equations  as  an 
initial  condition.  For  this  pod  a  sideslip  of  -0.5  degrees 
was  assumed.  The  resulting  lateral  motion  is  not  shown 
because  the  main  emphasis  is  on  the  longitudinal  motion.  The 
pod  is  unstable  laterally,  also,  but  the  initial  movement  is 
small  due  to  the  relatively  small  initial  conditions.  In 
addition  to  the  yaw  angle  and  rate,  if  the  ejection  rack 
provided  an  initial  roll  rate,  such  as  when  the  pod  is  loaded 
off-center  of  the  bomb  rack,  then  that  influence  could  also  be 
included. 

D.   PARAMETER  VARIATION 

The   parameter   variation   due   to   non-axisymmetry   is 
difficult  to  handle.    Comparison  with  shapes  with  known 
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(experimentally  obtained)   coefficients  might  provide  some 
degree  of  accuracy  in  predicting  more  accurate  results . 

E.   PRELIMINARY  OBSERVATION 

The  results  of  the  method  in  Section  II  were  a  vertical 
velocity  of  29.9  feet/sec  and  an  angular  velocity  of  64.1 
deg/sec.  The  aerodynamic  method  of  Section  III  begins  with  a 
vertical  velocity  of  29.9  fps,  quickly  diverging  to  an 
extreme  value.  The  pitch  rate,  q,  also  diverges  quickly. 
Because  the  pod  is  very  unstable  any  linear  approximation  of 
its  behavior  will  have  strict  limits .  The  coefficients  used 
in  this  calculation,  from  the  MISDATCOM  code,  are  not  valid 
beyond  a  fraction  of  the  trajectory.  The  valid  time  range  of 
prediction  is  presented  in  Section  V. 
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IV.    "NEAR"  PROGRAM 

A .   BACKGROUND 

A  computer  prediction  method  was  developed  by  Nielsen 
Engineering  and  Research,  Inc.,  (NEAR)  under  contract  to  the 
United  States  Air  Force.  The  work  was  performed  during  the 
period  1968  to  1972.  The  final  result  is  a  method  for 
predicting  the  six  degree-of-f reedom  store  separation 
trajectory  at  speeds  up  to  the  subsonic  critical  Mach  number. 
After  delivery  of  the  program  to  the  Government,  many  new 
capabilities  were  added.  The  code  used  for  this  paper  was 
obtained  from  the  Naval  Weapons  Center,  China  Lake,  Ca.  The 
program  has  been  widely  accepted  by  industry  and  government. 
The  code  encompasses  9, 900  lines  of  code  and  thus  requires  a 
device  with  sufficient  computer  memory  for  operation.  (This 
also  depends  upon  the  program  application.) 

The  aircraft  fuselage,  separated  store  body,  and  adjacent 
stores  are  modeled  using  point  sources  and  sinks.  Angle  of 
attack  effects  are  included  using  a  cross-flow  model .  The 
aircraft  wing  and  wing  pylons  are  modeled  using  planar  vortex 
lattice  models  which  include  dihedral,  camber,  and  twist  of 
the  aircraft  wing.  Thickness  strips  are  used  to  model  the 
thickness  of  the  aircraft  wing  and  pylons.  [Ref.  10 :p.  807] 
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The  capability  exists  to  install  multiple  sets  of  wings, 
fins,  or  canards,  and  to  use  active  control  surfaces  on  the 
store  by  inputting  the  control  laws  into  the  program.   Powered 
separations  may  be  simulated  by  inputting  the  store  thrust 
characteristics . 

The  NEAR  program  actually  consists  of  two  separate 
programs:  the  source  program  and  the  trajectory  program.  Both 
are  described  in  Section  IV. B  below. 

Alternate  separation  programs  include  USTORE  and  USSAERO 
codes.  USSAERO  was  developed  by  F .  A.  Woodward,  of  NASA,  as 
a  lower-order  panel  method.  USTORE  was  developed  from  USSAERO 
by  G.  J.  van  den  Broek,  of  the  National  Institute  for 
Aeronautics  and  Systems  Technology,  Pretoria,  South  Africa. 
[Ref.  ll:p.  309] 

B .   THEORY 

The  three  principal  tasks  in  the  prediction  of  a  store 
trajectory  are:  first,  the  determination  of  the  nonuniform 
flow  field  in  the  neighborhood  of  the  ejected  store;  second, 
the  determination  of  the  forces  and  moments  on  the  store  in 
this  flow  field;  and  third,  the  integration  of  the  equations 
of  motion  to  determine  the  store  trajectory  [Ref.  3]. 

1 .  Source  Program 

The  source  program  is  used  to  represent  an 
axisymmetric  body  as  a  distribution  of  sources  along  the  axis 
of  the  body.   It  provides  point  source-sink  distributions  to 
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represent  the  fuselage,  rack,  and  store  volumes.  The  program 
calculates  and  prints  the  source  strengths  and  locations. 
These  quantities  are  then  used  as  input  data  to  the  trajectory 
program. 

The  source  program  is  used  for  the  generation  of  an 
aircraft  or  pylon  model  which  is  then  used  for  the  trajectory 
program.  These  models  are  Mach  number  dependent  and  thus  need 
to  be  generated  for  each  test  case  with  different  Mach  number. 
For  ongoing  store  separation  studies,  a  good  database  of 
aircraft  models  is  required  and  should  be  available  from 
appropriate  Government  facilities.  The  F/A-18  model  was 
obtained  from  China  Lake  along  with  the  program  code.  This 
model  included  pylon  stations. 
2.  Trajectory  Program 

The  trajectory  program  uses  the  source-sink 
distributions  from  the  source  program,  and  additional 
information  to  first  determine  the  vorticity  distribution 
which  represents  the  wing-pylon  loading  including  interference 
of  the  fuselage,  rack,  and  stores,  and  then  to  calculate  the 
store  trajectory.  See  Sections  IV.C.3.d.(l)  and  (2)  for 
descriptions  of  the  vortex  lattice  method  and  the  panel 
method. 

Once  the  trajectory  program  input  has  been  generated 
for  a  single  flight  condition,  it  is  relatively  easy  for  the 
user  to  input  different  stores.   The  major  effort  is  to  obtain 
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the  initial  aircraft  fuselage,  wing,  and  pylon  models  for  this 
input . 

3 .  Shape  Modeling 
a .  Aircraft 

The  aircraft  fuselage  geometry  was  modeled  by  using 
a  Fortran  program  called  NGDELX  developed  by  L.  Gleason  of  the 
Naval  Weapons  Center.  This  program  provides  the  coefficients 
to  the  NEAR  polynomial  representation  of  equivalent  body.  The 
shape  is  divided  into  appropriate  segments.  Radius  values  for 
the  segment  points  and  maximum  radius  values  are  entered. 
These  coefficients  are  then  entered  into  the  source  program  to 
generate  the  source  representation  of  the  fuselage.  Figure  9 
shows  the  coefficients  for  the  F/A-18  fuselage.  Following  is 
the  equation  used  to  calculate  the  coefficients: 

r.q/L  =  Cj  +  C7*(C2*(X/L)2  +  C3*(X/L)+C4)  +  C5*  (X/L)  +  C6*  (X/L) 2 


Jb.  Pylons 

The  pylons  used  are  modeled  the  same  as  in  Section 
IV.C.l.a  above.  The  geometry  is  used  to  obtain  the  NEAR 
coefficients . 
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Polynomial  Rapraatntatlon  of  Equivalent  Body 
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Figure  9    F/A18  NEAR  Polynomials 

c.  Pod  or  Missile 

Again,  the  procedure  in  Section  IV.C.l.a  is  used  to 
represent  the  shape  as  a  set  of  polynomial  coefficients.  The 
NGDELX  program  is  again  used  to  represent  the  missile  shape  as 
an  "equivalent  body  of  revolution"  (EBR)  .  NGDELX  is  very  easy 
to  use,  but,  again,  only  represents  axisymmetric  shapes. 

d.  Wing    (Lifting  Sv~£e.ces) 

The  aircraft  wings,  missile  fins,  or  any  other 
lifting  surfaces  are  modeled  using  a  vortex  lattice  method  to 
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represent  the  loaded  wing.    The  thickness  is  also  modeled 
using  the  lattice  method. 

(1)  Vortex  Lattice  Methods.  There  are  several 
variations  of  the  vortex  lattice  method  that  are  presently 
available  and  have  proven  to  be  very  practical  and  versatile 
theoretical  tools  for  the  aerodynamic  analysis  and  design  of 
planar  and  non-planar  configurations.  [Ref.  12  p.  27]  The 
vortex  lattice  method  represents  the  wing  as  a  planar  surface 
on  which  a  grid  of  horseshoe  vortices  is  superimposed.  The 
velocities  induced  by  each  horseshoe  vortex  at  a  specified 
control  point  are  calculated  using  the  law  of  Biot-Savart .  A 
summation  is  performed  for  all  control  points  on  the  wing  to 
produce  a  set  of  linear  algebraic  equations  for  the  horseshoe 
vortex  strengths  that  satisfy  the  boundary  condition  of  no 
flow  through  the  wing.  The  vortex  strengths  are  related  to 
the  circulation  and  the  pressure  differential  between  the 
upper  and  lower  wing  surfaces.  The  pressure  differentials  are 
integrated  to  yield  the  total  forces  and  moments.  [Ref.  13 
p. 261]  For  a  rigorous  introduction  to  the  vortex  lattice 
method,  see  Reference  13. 

(2)  Panel  Methods.  The  configuration  is 
modeled  by  a  large  number  of  elementary  quadrilateral  panels 
lying  either  on  the  actual  aircraft  surface,  or  on  some  mean 
surface,  or  on  a  combination  thereof.  To  each  elementary 
panel,  there  is  attached  one  or  more  types  of  singularity 
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distributions,  such  as  sources  ,  vortices,  and  doublets. 
These  singularities  are  determined  by  specifying  some 
functional  variation  across  the  panel  (e.g.,  constant,  linear, 
quadratic,  etc.),  whose  actual  value  is  set  by  corresponding 
strength  parameters .  These  strength  parameters  are  determined 
by  solving  the  appropriate  boundary  condition  equations.  Once 
the  singularity  strengths  have  been  determined,  the  velocity 
field  and  the  pressure  field  can  be  computed.  [Ref .  13.*p.  258- 
259] 

C .   PROCEDURE 

1 .   Source  Program 
a .   Input 

The  program  input  consists  of  the  polynomial 
representation  of  the  equivalent  body  of  revolution  (EBR) 
values  obtained  in  Section  IV.B.3.a.  The  surface  of  the  EBR 
is  approximated  by  these  polynomials,  which  represent  the  EBR 
x,r  distribution.  The  source  program  is  then  run  with  a  user- 
specified  finite  distribution.  The  surface  obtained  from  this 
source  distribution  can  then  be  compared  with  the  polynomial 
surface.  Additional  runs  may  be  required  to  closely  match  the 
surfaces.  The  aircraft  fuselage,  pylon (s),  and  the  adjacent 
store  EBR  values  are  input  along  with  two  variables,  NRAT  and 
PERCR.  NRAT  is  the  number  of  segments  which  the  body  will  be 
divided  for  the  specification  of  the  source  distribution. 
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PERCR  is  the  source  spacing  for  each  NRAT  segment  and  is  input 
as  a  fraction  of  the  local  body  radius  of  the  segment . 

Adjacent  stores  are  included  in  the  aircraft  source 
input.  In  this  way  it  becomes  part  of  the  aircraft 
configuration.  The  separated  store  is  not  input  to  the  source 
program  since  the  trajectory  program  calls  the  source  program 
during  its  execution  to  model  the  store.  this  is  due  to  the 
fact  the  store  changes  position  during  the  program  run  and  the 
sources/sinks  will  change  values.  Figure  10  lists  a  sample 
format  for  the  source  program  input . 
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Figure  10    Source  Program  Input 
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b.     Output 

The  source  program  output  is  used  as  input  to  the 
trajectory  program.  However,  it  is  not  directly  read  into  the 
trajectory  input  and  must  be  entered  via  keyboard.  The  form 
of  the  output  is  shown  as  a  partial  output  in  Figure  11. 
Source  locations  and  strengths  are  listed. 
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Figure  11    Source  Program  Output 
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2  .   Trajectory  Program 

The  trajectory  program  uses  the  output  of  the  source 
program,  and  the  input  of  the  wing  configuration  and  pylons, 
along  with  store  information  to  calculate  the  trajectory  of 
the  separated  store.  It  is  not  the  intent  of  this  report  to 
fully  explain  the  aircraft  modeling  details  of  the  program. 
The  assumption  is  made  that  the  appropriate  aircraft  model  is 
available  for  the  correct  flight  condition  and  configuration. 
The  applicable  store  inputs  to  the  program  are  discussed  and 
the  resulting  trajectory  examined.  Nielson  contains  a 
thorough  description  of  the  theory  and  input  of  the  aircraft 
model  in  Reference  3 . 

a .  Input 

The  example  input  to  the  trajectory  program  was 
1064  lines  in  size.  The  store  data  entry  begins  at 
approximately  980  lines,  depending  upon  aircraft  input.  The 
store  input  will  vary  according  to  the  physical  configuration 
of  the  store.  For  the  example  separation,  Input  Items  35-43 
and  48-51  were  used  to  input  store  physical  data.  The  ejector 
information  was  input  in  Items  44A-47.  This  is  where  the 
fifth-order  polynomial  representation  of  the  ejector  force  was 
used.  Data  concerning  initial  and  final  times,  and 
integration  step  sizes  were  entered  in  Input  Item  72 .  The 
correct  input  form  of  the  above  data  is  contained  in  the  NEAR 
code  users  manual  [Ref.  14]  and  is  not  presented  here. 
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b .  Output 

The  output  of  the  trajectory  program  is  dependent 
upon  the  input  conditions,  i.e.,  the  number  of  stores,  number 
of  store  segments,  and  number  of  integration  steps,  etc.  The 
example  case  output  is  4,524  lines.  The  first  1,840  lines  of 
output  are  aircraft-specific  output,  such  as  the  source  and 
vorticity  information. 

The  last  2,650  lines  of  code  contain  the  trajectory 
data.  The  ejector  force  data  is  repeated  as  is  the  store 
shape  data,  reference  dimensions  and  inertial  characteristics. 

Proceeding  the  ejector  data,  a  data  block  is 
presented  for  each  time  integration  step.  This  data  block 
contains  the  parameters  which  describe  the  store  displacement, 
velocity,  and  acceleration  during  the  separation.  A  partial 
sample  output  of  this  section  of  the  output  file  is  included 
in  Figure  12.  Plotting  routines  can  be  developed  to 
facilitate  data  analysis  of  the  simulation  but  have  not  been 
developed  for  this  study.  The  following  are  short 
descriptions  of  some  of  these  parameters .  A  comparison  with 
the  results  of  the  method  of  Section  III  follows  in  Section 
IV.  D. 

The  force  and  moment  coefficients  are  listed  by 
their  individual  contributions.  The  effects  of  buoyancy, 
slender  body  theory,  crossflow,  and  empennage  are  totaled  to 
give  an  effective  CN,  CY,  C^,  and  Cln. 
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Figure    13         Trajectory   Output     (Cont) 
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The  ejector  force  and  moment  components  are  listed. 
The  vertical  force  history  matches  the  force  profile  from 
Section  II  very  closely.  The  store's  nose,  inertial  reference 
center,  and  the  base  position  in  the  fuselage  coordinate 
system  is  listed.  The  separation  distance  from  the  initial 
position  is  also  listed  as  output .  Plots  of  these  parameters 
will  indicate  the  miss  distance  of  the  shape  from  the  wing. 
The  store  translational  and  rotational  velocities  are  output 
in  both  the  store's  reference  frame  and  the  fuselage  reference 
frame . 
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V.     RESULTS  AND  COMPARISONS 

A.  EJECTOR  FORCES  AND  MOMENTS 
1 .  Vertical  Velocity 

The  predicted  end-of-stroke  vertical  velocity  from  the 
hand-calculation  was  29.9  fps,  comparing  with  the  NEAR  code 
results  of  30.84.  This  is  a  3%  difference.  The  predicted 
value  from  the  ejector  graph  of  Reference  4  yields  24.5  fps. 
The  ejector  graph  is  representative  of  the  average  store  used 
in  the  Naval  inventory,  which  more  than  likely  would  have 
control  surfaces  for  stability  reasons.  These  surfaces  would 
also  provide  aerodynamic  damping  which  could  be  the  source  of 
the  discrepancy. 

2  .  Pitch  Rate  Q 

The  predicted  pitch  rate,  q,  of  1.12  radians/second 
matches  well  with  the  NEAR  end-of-stroke  q  of  1.08 
radians/second.  For  most  applications,  this  correlation  of  4% 
should  have  sufficient  accuracy. 

B.  TRAJECTORY  COMPARISON 

The  main  consideration  in  determining  the  accuracy  of 
these  trajectory  methods  is  the  safe  separation  (or  jettison) 
of  the  pod.    Vertical  distance  and  pitch  attitude  are 
indications  of  the  pod  separation  while  velocity  and  angular 
velocity  are  indications  of  the  store  loads.    Therefore, 
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criteria  for  judging  the  correlation  between  the  linear  method 
and  the  NEAR  code  must  be  based  on  the  practical  standpoint  of 
a  reasonable  assumption  of  safe  separation.  The  individual 
parameters  are  addressed  below.  For  determining  the  valid 
time  range  of  the  aerodynamic  prediction  method,  the  NEAR  code 
was  used  as  reference . 

1 .  Vertical  Separation  Z 

Vertical  distance  z  in  the  linear  method  is  derived 
from  the  vertical  velocity  term.  Therefore  the  error 
(deviation  from  NEAR  results)  will  be  slightly  less  due  to  the 
integration  effects.  However,  given  the  initial  vertical 
velocity,  the  NEAR  code  predicts  a  five  foot  separation.  The 
error  limit  must  be  set  as  an  acceptable  variation  of  this 
distance.  Based  upon  practical  knowledge  of  store  separation, 
a  20%  error  limit  was  set.  Figure  14  compares  the  linear 
aerodynamic  method  results  comparison  with  the  NEAR  results  in 
the  amplified  time  region  of  the  first  0.2  seconds. 

2  .  Pitch  Angle 

Since  the  pitch  angle  is  used  solely  for  determining 
the  safe  separation  and  is  not  used  for  more  accurate 
applications,  like  rocket  engine  ignition,  the  effect  of  an 
error  in  pitch  will  vary  as  the  sine  of  the  angle  deviation. 
This  could  result  in  an  error,  due  to  rotation  about  the 
center  of  gravity,  of  an  order  of  7.2  feet  times  the  sine  of 
the  error  angle.    The  pitch  angle  error  must  be  fairly  small 
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initially,  but  can  grow  as  time  increases  due  to  increased 
separation  distance  and  its  decreasing  effect.  Therefore,  a 
conservative  error  limit  is  set  for  0.1  sec.  At  0.1  seconds, 
the  NEAR  pod  has  rotated  4.5  degrees.  A  30%  error  would  give 
a  vertical  error  of  0.2  feet,  which  is  approximately  10%  of 
the  actual  separation.  Figure  15  shows  a  comparison  of  the 
linear  method  with  the  NEAR  results.  Figure  16  provides  the 
error  between  the  methods  as  a  function  of  time. 
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Figure  16    Parameter  Error 
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3 .  Pitch  Rate 

The  required  accuracy  for  pitch  rate  will  be 
determined  by  its  application.  For  instance,  navigational 
considerations  will  require  quite  accurate  predictions  of  the 
angular  loading  of  its  components.  However,  for  this  case, 
there  is  not  a  requirement  of  this  kind.  The  divergent  nature 
of  the  pitch  rate  from  the  linear  aerodynamic  model  does  not 
present  a  problem  for  this  case,  beyond  the  fact  that  the 
pitch  angle  follows  this  trend.  Therefore,  an  error  limit  is 
not  placed  upon  the  pitch  rate. 

4 .  Vortical  Velocity 

Vertical  velocity  is  divergent  and  represents  the  lack 
of  aerodynamic  damping  in  the  linear  model.  Therefore,  it 
does  not  provide  meaningful  information  beyond  the  phase 
relationship  with  the  vertical  separation  distance.  For  this 
reason,  an  error  limit  is  not  placed  upon  the  velocity. 

C.   AUGMENTED  AEROPREDICTION 

Due  to  the  limited  valid  time  limits  on  the  aerodynamic 
method,  an  alternative  approach  was  investigated.  Obviously, 
the  non-linear  operating  region  is  reached  very  quickly.  The 
two  dominant  terms  of  the  A-matrix  are  the  damping  terms,  C^ 
and  C^.  The  unstable  pod  is  shown  from  the  NEAR  code  to 
quickly  reach  a  somewhat  steady-state  value  of  pitch  rate  as 
it  tumbles.  An  effort  was  made  to  calculate  a  new  (and  valid) 
value  of  Cng.   This  approach  modeled  the  tumbling  store  as  a 
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cylinder  of  constant  radius.   Using  a  value  of  Cd=1.0,  and  the 
end-of-stroke  pitch  rate  of  1.12  radians/second,  new  values  of 
damping  coefficients  were  calculated. 
1  .  C^-  Coefficient 

The  new  pitch  damping  coefficient  of  -.1054  was  used 
versus  the  MISDATCOM  value  of  -.04.  This  is  an  obvious 
improvement  and  also  realistic  since  the  pod  will  definitely 
experience  aerodynamic  damping,  due  to  pitch  rate,  that  would 
not  have  been  predicted  by  MISDATCOM. 

2 .  C^-  Coefficient 

The  value  of  the  aerodynamic  damping  effects  of  change 
in  angle  of  attack,  C^,  was  combined  with  C^  in  the  MISDATCOM 
output.  Therefore,  the  valid  value  of  C^  is  unclear.  For 
this  model  of  a  rotating  pod,  it  was  assumed  that  C^  would 
have  minimal  effect  on  the  damping  and  a  value  of  0.0  was 
assigned. 

3 .  Trajectory  Prediction 

The  predicted  parameters  using  the  above  coefficients 
are  contained  in  Figures  17,  18  and  19.  As  shown  in  the 
figures,  there  is  considerable  improvement  in  the  correlation 
with  the  NEAR  trajectories  for  pitch  angle  and  pitch  rate. 
The  more  realistic  value  for  C^  obviously  improves  the 
results.  The  time  range  of  validity  for  these  two  parameters 
is  now  estimated  to  be  1.0  seconds.  However,  the  most 
important  parameter,  vertical  separation  distance,  is  still 
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outside  of  reasonable  accuracy  limits.  Therefore,  the 
augmented  plant  matrix  method  is  still  unusable  for  complete 
trajectory  prediction.  Figure  20  shows  the  parameter  error 
for  the  augmented  system. 

Although  there  is  substantial  correspondence  in  the 
partial  comparison  between  the  NEAR  results  and  these  results, 
there  are  not  enough  data  points  to  substantiate  this  method. 
Further  work  using  different  shapes  need  to  be  done. 
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Figure  17    Vertical  Separation  Distance  (Augmented) 
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Figure  18    Pitch  Angle  (Augmented) 
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Figure  19    Pitch  Rate  (Augmented) 
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Figure  20    Parameter  Error  (Augmented) 
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VI  .   CONCLUSIONS 

In  the  preliminary  design  application,  MISDATCOM  provides 
an  easy-to-use  method  to  predict  moderately  accurate  store 
aerodynamic  characteristics.  These  aerodynamic  coefficients 
can  be  used  to  model  the  store  dynamics  for  design  work. 
However,  their  use  in  the  linearized  equations  of  motion  for 
an  externally  ejected  store  trajectory  prediction  is  highly 
limited. 

Comparisons  with  the  Nielson  Engineering  and  Research 
(NEAR)  code  show  that  trajectory  values  are  outside  of 
reasonable  accuracy  limits  for  use  in  store  separations.  The 
store  separation  vertical  distance  and  pitch  angle  values 
diverged  within  0.1  seconds.  Therefore,  the  use  of  the 
linearized  equations  of  motion  cannot  provide  useful 
trajectory  information. 

However,  a  modification  to  these  coefficients  to  more 
accurately  reflect  the  existing  aerodynamic  damping  does 
provide  correlations  with  the  NEAR  code  of  within  30%  for 
pitch  angle  and  10%  for  pitch  rate.  Use  of  these  "augmented" 
linear  equations  of  motion  for  attitude  prediction  shows 
sufficient  accuracy  for  many  applications. 

The  NEAR  code  is  widely  used  for  ejected  store  trajectory 
prediction  in  the  aviation  field.   Its  use  for  safe  separation 
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investigation  should  be  considered  as  standard  for  obtaining 
accurate  results  .  The  cost  of  the  Nielson  code  in  both 
engineering  time  and  computational  time  increases  with  the 
level  of  information  produced,  but  the  need  for  safe 
separation  accuracy,  aerodynamic  load  distributions,  and 
vortex-induced  effects  make  the  extra  effort  justified. 

For  limited  applications  requiring  less  accuracy,  the 
augmented  linear  equations  of  motion  can  provide  some 
information  for  the  initial  trajectory  movement.  A  potential 
application  would  be  the  investigation  of  the  effects  of  small 
weight  or  inertial  modifications  to  existing  stores . 
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APPENDIX  A 

A.     POD  CHARACTERISTICS 

The  pod  in  question  was  an  ALE-29B  pod  modified  to 
accommodate  the  SUU-53  chaff  dispenser  and  associated 
equipment.  The  ALE-2  9  pod  was  flight  approved  by  the  Aircraft 
Configuration  Control  Board  (ACCB)  document  #75-004.  It  was 
originally  an  ALE-2  pod.  No  significant  aerodynamic  changes 
have  been  made  to  the  pod. 

The  pod  physical  characteristics  are  listed  in  Table  1. 
TABLE  1     POD  CHARACTERISTICS 


Length 

153. 

8  inches 

Radius    -minimum 

9.75 

inches 

-maximum 

11.20 

inches 

CG  location  -empty 

62.4 

inches  from  nose 

-full 

65.4 

inches 

Weight    -empty 

326 

lbs 

-full 

371 

lbs 

Moment  of  Inertias 

-Ixx 

300 

lb-ft2 

-Iyy 

1675 

lb-ft2 

-Izz 

1675 

lb-ft2 
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APPENDIX  B 


A.     FORCE  EQUATIONS 


IF  =  M*V  -  d/dt (M*V)  =   G 


Separating  variables 


IFX    =    Gx  lFy    =    Gy  IFZ    =    Gz 


Assuming    £Fxy    =    0, 

Integrating    the    remaining   equation, 

JlFz    dt    =    G2    -    Gj    =     (M*VZ)2    -     (M*V2)i 
Assuming   straight    and   level    flight,    Vzl   =    0    yields: 

JlFz   dt    =   M*AVZ    =   M*VZ 

Integrating  the  ejector  force  polynomial  yields  a  total 
impulse  of  345  lbf-secs.   Substituting  (wt=371  lb), 

345  lbf.secs  =  371  lbm  *  Vz 

Vz  =  1.1773  (lbf-sec/lbm)  *  32.174  (lbm-f t/lbf-sec2) 

Vz  =  2  9.9  fps 
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B.  MOMENT    EQUATIONS 

IM0    =    r    x    IF    =    r    x    M*V    =    d/dt(r    x    M*V)     =    H0 

where  H0  is  the  angular  momentum  about  pt  0.   Integrating  the 
moments, 

JlM0  dt  =  H  +  r  x  H 

Separating  the  component  equations  yields, 

£MX  =  Ixx*d)x  -  (Iyy  -  Izz)*0)y*COz 
ZMy  =  Iyy*^  -  (I„  -  Ixx)  *G)Z*CDX 
IMZ  =  Izz*d)z  -  (Ixx  -  Iyy)  *caK*COy 

Assuming  the  reference  frame  coincides  with  the  principal 
axes,  Ixy,xz,yz  =  0,  and  assuming  coxz  =  0,  yields  the  following: 

JlMy  dt  =  Iyy  *  (by 

toy   =1.12  rad/sec 
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APPENDIX  C 

A.     LONGITUDINAL  EQUATIONS  OF  MOTION 

The  equations  of  motion  for  the  pod  can  be  derived  from 
Newton' s  Second  Law  of  motion,  which  states  that  the  summation 
of  all  external  forces  acting  on  a  body  must  be  equal  to  the 
time  rate  of  change  of  the  momentum  of  the  body,  and  the 
summation  of  the  external  moments  acting  on  a  body  must  be 
equal  to  the  time  rate  of  change  of  the  moment  of  momentum 
(angular  momentum) .  The  time  rates  of  change  are  all  taken 
with  respect  to  body  coordinates  space.  [Ref.8] 


w  ,c%1 


y..v.-cY 


Figure  21    Store  Reference  Frame 
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£  AFx=/7?*  (U+W*Q-V*R) 

Y^  AFz=/7?*  (W+V*P-U*Q) 

£  &M=Q*Iy  +  P*R*  ilx"lg)  +  iP'-R2)  *JX 


£  AFx=m*  (C/+^*(?-V*i2) 

£  AF2  =  m*  (^+V*P-(7*0) 

£AM=£)*Jy 


By  restricting  the  disturbances  to  small  perturbations 
about  the  equilibrium  condition,  the  product  of  the  variations 
will  be  small  in  comparison  with  the  variations  and  can  be 
neglected,  and  the  small  angle  assumptions  can  be  made 
relative  to  the  angles  between  the  equilibrium  and  disturbed 
axes  . 


J^  AFx=/n*ti 

J^  AFz=m*  (w-U*q)  =m*  (w-U*&) 

£  AM=Iy*g=Iy*fc 


Expanding  the  applied  forces  and  moments  in  terms  of  the 
changes  in  the  aerodynamic  and  gravitation  forces  and  moments 
using  the  total  differential  form  yields: 
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£  dFx=  (dFjdu)  *dlT+  (dFx/dw)  *dw+  (dFx/dW)  *dw+  (dFx/dG)  *dS-  {dFx/db[ 

£  dFz=  idFjdu)  *du+  (dFjdW)  *dfV-  {dFjdw)  *dw+  (dFz/d&)  *dO~  (dFz/db[ 

Y,  AM=  idM/dU)  *dU+  (dM/dW)  *dW+  {dM/dW)  *dW+  {dM/d&)  *d& 

Non-dimensionalizing  the  terms  in  the  above  equations 
yield  the  aerodynamic  coefficients.  Following  are  the 
definitions  of  the  longitudinal  stability  coefficients  and 
derivatives : 

Cxu      Variation  of  drag  and  thrust  with  u. 

Cxa     Lift  and  drag  variations  along  the  x-axis . 

Cw      Gravity 

Cx(i     Downwash  lag  on  drag. 

Cxq     Effect  of  pitch  rate  on  drag. 

Czu     Variation  of  normal  force  with  u. 

Cxa  Slope  of  the  normal  force  curve. 

Cxi     Downwash  lag  on  lift  of  tail. 

Czq     Effect  of  pitch  rate  on  lift. 

Cmu     Effects  of  thrust,  slipstream,  and  flexibility. 

Cma     Static  longitudinal  stability. 

Cm<i     Downwash  lag  on  moment . 

Cmq     Damping  in  pitch. 

Arranging  the  linearized  longitudinal  equations  into 
matrix  form  yields  the  following  matrix  equation.   The  control 
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matrix  B  is  not  shown  here  but  is  included  in  the  Control-C 
program.  Since  there  are  no  active  control  surfaces,  the  B- 
raatrix  is  a  zero  matrix. 


Au 

xu 

x. 

0 

-g 

Au 

Aw 

Zu 

Zw 

u 

0 

Aw 

Aq 

= 

MU+M„*ZU 

MW+MW*ZW 

Mq+M„*U 

0 

Aq 

AO. 

0 

0 

0 

0 

Ae 
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B.    LATERAL  EQUATIONS  OF  MOTION 

In  the  same  manner  used  to  obtain  the  longitudinal 
equations,  the  lateral  equations  are  derived.  The  same 
assumptions  using  perturbation  theory  apply.  Following  are 
the  equations  in  their  developing  order. 

£  AFy=/77*  (  V+U*R-W*P) 
£  AL=P*IX-R*JXZ+Q*R*  UM~Iy)  ~P*Q*JXZ 
£  AN=R*IZ+P*JXZ+P*Q*  (Iy-Ix)  +  Q*R*JXZ 

^T  Afr=m*  ( v+U*r+u*r) 
Y,*N=r*Iz-p*Jxz 

^2  AFy=/77*  ( v+U*r) 
J2&L=p*Ix-r*Jxz 
Y;*N=r*I2-p*Jxz 

£  dFy=  {dFy/dB)  *dB+  (dFy/dlf)  *dT+OFy/3*)  *d«+  (8Fy/6$)  *d*-  (dFy/- 


Ap 

Yp/u 

Yp/u 

-(l-Yr/u) 

g*cos0/u 

Af 

Ap 

Lp 

Lp 

Lr 

0 

Ap 

Ar 

= 

Np 

Np 

Nr 

0 

Ar 

A<t> 

L0 

1 

0 

0 

A(D. 
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APPENDIX    D 


m- Input (' ente r    m    ( lbs ) : ' )/3 2 . 2 ; 

Uo-lnput ( ' ente r    Uo    (fpB)t')r 

S-input( 'enter    S    ( i n* *  2 ) t * ) /l 4 4 ; 

b-inpuM 'entet    b    (ln);')/12; 

Ro-lnput( 'enter  Ro  ( slugs/ft ** 3  )  t ')  ; 

cref-lnput (' enter  c  rtf I in ) i ' J/12 ; 

Ix-input( 'enter  Ixx  ( ft-lb) j ' )/32. 2| 

Iy-input( 'enter  IYY  ( f t-lb ) i ' )/32 . 2 ; 

Iz-lnput( 'enter  Izi  ( ft-lb) t • )/32 . 2 | 

g-32.174; 

Cn-lnput (' enter  Cn  (Normal  force  coefficient)!'); 

Cm-1 nput (' enter  Cm  (Pitching  Moment  Coeff lclent )  I  '  )  | 

Ca-lnput (' enter  Ca  (Axial  Force  Coefficient)!')) 

Cy-tnput (' enter  Cy  (Side  rorce  Coefficient ) t ' ) ; 

CLN-lnput( ' enter  CLN  (Yawing  Moment  Coeff):'); 

CLL-lnput( ' enter  CLL  (Boiling  Moment  Coeff);'); 

CNA-lnput (' enter  CNA  (Normal  force  Derlv  wet  AOA):'); 

CMA-lnput (' enter  CMA  (Pitching  Moment  Derlv  wrt  AOA ) t ' ) ; 

CYB-input( ' enter  CYB  (Side  Force  Derlv  wrt  Beta);'); 

CLNB-lnput( ' enter  CLN8  (yawing  Moment  derlv  wrt  Beta);'); 

CLLB-lnput( ' enter  CLLB  (Rolling  Moment  deriv  wrt  Beta);'); 

CL-input( 'enter  CI  (lift  coeff);'); 

CD-Input (' enter  Cd  (Drag  coeff);'); 

CNQ-lnput( 'enter  CNQ  (Normal  Force  Deriv  wrt  Q)  ;  '  )  ; 

CMQ-lnput( 'enter  CMQ  (Pitching  Moment  Derlv  wrt  0);'); 

CNAOADOT-lnput( 'enter  CN  AOAdot  (Normal  Force  Derlv  wrt  AOAdot);'); 

CMAOADOT-lnput( 'enter  CM  AOAdot  (Pitching  Mom  Deriv  wrt  AOAdot);'); 

q- . 5*  ro*uo*uo; 

Xu--2*Cd*q*s/(m*uo) ; 

Xw-Cl*q*s/(m*uo) ; 

Zu — 2*Cl*q*s/(m*uo) ; 

Zw--Cd*q*s/(m*uo ) ; 

Zwd-Cnaoadot*cref*a*q/(2*uo»uo*m); 

Mw-Cma*q*B*cref/( Iy*uo) ; 

Mu-0; 

Ma-uo*Mw; 

Mwd-Cmaoadot *cref *q*a*cre  f/(2*uo*uo*Iy); 

Mad-uo*Mwd; 

Mq-Cmq*cref *q*S*cref/( 2*uo* Iy ) ; 

A-(Xu  Xw  0  -g;Zu  Zw  uo  0; ( Mu+Mwd*Zu )/( 1-Zwd )  ( Mw*Mwd*Zw)/( 1-Zwd) . . . 

(Mq+Mwd*uo)/( 1-Zwd)  0;0  0  1  0); 
B-(0  0  0  OjO  0  0  0J';C-10  1  0  0);D-[0  0); 
1C-10  0  0  0)'  ; 

IC(2)-input( 'enter  initial  vertical  velocity  (fps);'); 
IC( 3)-input< 'enter  initial  pitching  rate  (rad/s);'); 
ttt-input( ' enter  end  time  (sec)'); 
t-0;  .001ittt| 
uuu-(ttt/.001)+l> 
u-0'ones ( 2 , uuu ) ; 
simu(  '  IC  ,  IC)  ; 
y-slmu(a,brc,d,u,t) ; 
erase ; 
plot(t.y); 


Figure  22    Control-C  Program  For  EOM  Simulation 
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Figure    23         MISDATCOM    INPUT/OUTPUT 
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Figure    24         MISDATCOM    Input/Output     (Cont) 
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Figure    25         MISDATCOM    Input /Output     (Cont) 
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